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1 Introduction 

The study of critical phenomena and phase transitions is a traditional subject 
of statistical physics which has known its "modern age" , since powerful ap- 
proaches have been developed (renormalization group, conformal invariance, 
sophisticated simulation algorithms, . . . ). Simplified models attracted a lot 
of attention. This is essentially due to a spectacular property of continuous 
phase transitions at their critical point, scale invariance, which leads to an 
extreme robustness of some quantities, like the critical exponents which are 
thus referred to as universal quantities. Only very general properties (space 
or spin dimension, symmetry, range of interaction, . . . ) determine the uni- 
versality class. This makes the theory of critical phenomena a very efficient 
and predictive tool: As soon as one knows the general characteristics of a 
physical system from general symmetry arguments, it is possible in principle 
to predict exactly the "shape" of the singularities which are developed at the 
critical point. The term "exact" is here understood rigorously, for example 
a two-dimensional system with the symmetries of an Ising model, should it 
be a magnet, an alloy or anything else, will exhibit a diverging susceptibility 
X ~ |r — Tcl"^/** (see a sketch in Fig. [1]) with the precise value 7/4 for the 
exponent. 

On a theoretical ground, the major two-dimensional problems (Ising 
model and its generalizations like the Potts model and the percolation model. 
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Fig. 1. Typical behaviour of the susceptibihty at a second order phase transition. 
The quantities 7 = 7' and 1+ /F- are universal. 

XY model, Heisenberg model and so on) are essentially solved at least for 
their critical singularities (when a second-order phase transition is indeed 
present), but critical exponents are not the only universal quantities at a 
critical point. The universal character of appropriate combinations of crit- 
ical amplitudes [1] is also an important prediction of scaling theory but in 
some cases these combinations remain uncompletely determined and subject 
to controversies. 

The Potts model [2,3], as one of the paradigmatic models exhibiting 
continuous phase transitions is a good frame to consider the question of 
universal combinations of amplitudes. The universality class of the Potts 
model at its critical point is parametrized by the number of states q. The 
two-dimensional Potts model with three and four states can be experimen- 
tally realized as strongly chemisorbed atomic adsorbates on metallic sur- 
faces at sub- monolayer concentrations [4]. Although critical exponents could 
be measured quite accurately for adsorbed sub-monolayers, confirming that 
these systems actually belong to the three-state [5] or to the four-state Potts 
model classes [6], it is unlikely that the low temperature LEED results can 
be pushed [7] to determine also the critical amplitudes. Therefore, the nu- 
merical analysis of these models is the only available tool to check analytic 
predictions. 

The critical amplitudes and critical exponents describe the behaviour of 
the magnetization m, the susceptibility x^ the specific heat C and the corre- 
lation length ^ for a spin system in zero external fielcQ in the vicinity of the 
critical point 

^ In this paper we only deal with the physical properties in zero magnetic field. 
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M(t) W B_{-Tf, T < 0, (1) 

X(r)«r±|r|-^ (2) 
Xt(t) « rT{-Tr\ r < 0, (3) 

dr) « ^kl-, (4) 

^W^^o^kr^ (5) 

Here r is the reduced teraperature t — {T — Tc) /T and the labels ± refer to 
the high-temperature and low-temperature sides of the critical temperature 
Tc- For the Potts models with g > 2 a transverse susceptibility xt can be 
defined in the low-temperature phas^. 

Critical exponents are known exactly for 2D Potts model [8-11] through 
the relation Xe = {1 — oi)lv to the thermal scaling dimension 



1+2/ 

and the relation Xfj = (3/v to the magnetic scaling dimension 



(6) 



(7) 



where the parameter y is related to the number of states q of the Potts 
variable by the expression 

cos^ = i^/g (8) 

The central charge of the corresponding conformal field theory is also simply 
expressed [9] in terms of y 

0^1 -i^l. (9) 

Analytical estimates of critical amplitude ratios for the q-state Potts mod- 
els with g = 1, 2, 3, and 4 were recently obtained by Delfino and Cardy [12]. 
They used the two-dimensional scattering field theory of Chim and Zamolod- 
chikov [13] and estimated the central charge c = 0.985 for 4-state Potts model, 
for which the exactly known value is c = 1. Reporting these approximate val- 
ues in ([9]), one can calculate the scaling dimensions from ©-((T]) and get the 
values Xa = 0.13016 and x^ = 0.577, to be compared respectively to the exact 
values 1/8 and 1/2. The discrepancy is around 4 and 15 per cent, emphasiz- 
ing the difficulty of the 5 = 4 case (to give an idea, in the case of the 3-state 
Potts model, a similar analysis leads to a very good agreement with less than 
one percent deviation). 

The universal susceptibility amplitude ratios F^/ and Ft I F^ were also 
calculated in [12] and [14]. The figures obtained are the following. 



* In the following we will use equally the notations Fl or _r_ for the longitudinal 
susceptibility amplitude in the low temperature phase. 
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q = 3: F+IFl = 13.848, rr/A, = 0.327, 
g = 4 : r+/rL = 4.013, Tr/ri, = 0.129. 



(10) 
(11) 



These results have been confirmed numerically in the case g = 3 by several 
groups, r+/rL ~ 10 and rr/TL « 0.333(7) in Ref. [14] (Monte Carlo (MC) 
simulations), F^/Fl = 14±1 in Ref. [15] (MC and series expansion (SE) data) 
and quite recently, these results were confirmed and substantially improved, 
F+IFl = 13.83(9), Ft/Fl = 0.325(2), by Enting and Guttmann [16] who 
analysed new longer series expansions. 

The 4-state Potts model was also studied through MC simulations in 
Ref. [14], but the authors considered that their data were not conclusive. An- 
other MC contribution is reported by Caselle, et al [17], F+/Fl = 3.14(70), 
and Enting and Guttmann [16] also analysed SE data for the 4-state Potts 
model and found F^/F^ = 3.5(4), Ft/Fj^ = 0.11(4) in relatively good agree- 
ment with the predictions of [12] and [14]. The situation thus seems to be 
clear, although the use of the logarithmic corrections in the fitting procedure 
of MC data was questioned, e.g. in [16]: [Caselle et al] estimat,es depend crit- 
ically on the assumed form of the sub-dominant terms, and on the further 
assumption that the other sub-dominant terms, which include powers of loga- 
rithms, powers of logarithms of logarithms etc, can all be neglected. We doubt 
that this is true. 

Let us recall that the existence of logarithmic corrections to scaling in 
the 4-state Potts model was pointed out in the pioneering works of Cardy, 
Nauenberg and Scalapino [18,19], where a set of non- linear RG equations 
were proposed. Their discussion was later extended by Salas and Sokal [20] . 

Generically, the logarithmic corrections appear as corrections to scaling. 
We mentioned above that in the vicinity of a critical point, a susceptibility 
for example diverges like %(''") ~ r±\T\~'^ . This is true, but this singular 
behaviour can be superimposed to a regular signal (e.g. Da + Di\t\ + . . .), and 
the leading singular behaviour itself needs to be corrected when we consider 
the physical quantity away from the transition temperature. The expression 
for x(r) then takes a form which can become "terrific" : 

X(r) = Do -|- Dijrl + . . . regular backgromid 

+ i~'|T|~'''(l+ leading singularity 



Together with the amplitude and exponent associated to the leading singu- 
larity, F and 7, appear corrections to scaling due to the presence of irrelevant 
scaling fields (a^") and A, a'^"^ and A', . . . ), analytic corrections due to non- 
linearities of the relevant scaling fields (6^"^), or multiplicative logarithmic 



W^)|T|^+a(2)|T|2^-f .. 

+6W|t|+6(2)|^|2 + ...) 



leading corrections 



next corrections 



analytic corrections 
logarithmic corrections 
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corrections (<i and . . . ) . These logarithmic coefficients may have different 
origins (see e.g. in Ref. [1] and references therein). They can be due to the 
upper critical dimension, to poles in the expansion of regular and singular 
amplitudes, or to the presence of marginal scaling fields. The 4-state Potts 
model belongs to this latter category. 

Some of the quantities indicated above are universal. This is the case of 
the exponents as well as of many combinations of coefficients. In the present 
paper we are interested in the amplitude of the leading singular term, but its 
precise determination can be affected by the form of the corrections. We shall 
be concerned with the following universal combinations of critical amplitudes 

^ £± £l R+-^ (12) 

We present, for the 4-state Potts model, more accurate Monte Carlo data 
supplemented by a reanalysis of the extended series made available by Enting 
and Guttmann [16] and we address the following question: Is it possible to 
devise some procedure in which the role of these logarithmic corrections is 
properly taken into account? 



2 Amplitudes and universal combinations 

The scaling hypothesis states that the singular part of the free energy density 
can be written in terms of the deviation from the critical point, t = [T—Tc) /T 
and h = H — He, 

/si„g(r, h) = b-'^F^iKrby^T, Khby^h) (13) 

where F±{x, y) is a universal function (actually there is one universal function 
for each side r > or t < of the critical point) and Kr and ku are "metric 
factors" which contain all the non universal aspects of the critical behaviour. 
D is the space dimension. Let us stress that the functions F± are universal in 
the sense that some details of the model are irrelevant (e.g. the coordination 
number of the lattice (so long as it remains finite) , the presence of next nearest 
neighbour interactions, etc) but they depend on the boundary conditions or 
the shape of the system. The metric factors on the other hand depend on 
these details, and the universal combinations are obtained when the metric 
factors are eliminated from some combinations. 

The connection with scaling relations can be shown with an exam- 
ple. From Eq. (|13|). we also deduce similar homogeneous expressions for 
the magnetization, M(t, ft,) = b^-^^y^ KhM±{x,y) and the susceptibility, 
xiT,h) = b-'^+'^y'^KlX±{x,y). The choice b = (k^|t|)-i/?^- and ft = leads 
(for example below the transition temperature) for the following combination 
of quantities 
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CiT,0)x{r,0) 2_ ,.2_„_2;3_^ C_(1,0)X_(1,0) _ 
" m2(r,0) =("^""") " Af^(l,0) = ^'^^ 

The prefactor takes the value 1 thanks to the weU known scaHng rela- 
tion between critical exponents a + 2/3 + 7 = 2. Thus it follows that the 
above combination is a universal number. From the definition of magnetiza- 
tion, specific heat and susceptibility amplitudes in zero magnetic field, e.g. 

M(t,0) = \T\^^KhM_{l,0) = B^lrf by virtue of Eq. HI), this universal 
number is in fact a combination of amplitudes, R = A-F- / B^. We have 
similar universal combinations above the critical temperature or associated 
to other scaling relations. 



3 RG approach for the Potts model and logarithmic 
corrections at qr = 4 

Let us remind that the q-state Potts model is an extension of the usual lattice 
Ising model in which the site variables Si (abusively called spins) can have 
q different values, Si = 0, 1, . . . q — 1 but the nearest neighbour interaction 
energy —JSg^^sj only takes two possible values, e.g. —J and depending 
whether the neighbouring spins are in the same state or not. The Hamiltonian 
of the model reads as 

H = -jJ2Ss.s,- (15) 

(u> 

At the early times of real-space renormalization, the application to the pure 
Potts model led to some difficulties: the impossibility to affect a particular 
value for the spin of a cell after decimation due to a too large number of 
states (see Fig. [5]). 
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Fig. 2. Decimation of spin blocks for the Ising model (left) and high-g Potts model 
(right). In the latter case, the state of many cells cannot be decided by a simple 
majority rule. 



It was thus necessary to extend the parameter space, introducing new 
variables, called vacancies, to replace the question marks in Fig. ([2]) and 
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to study an annealed disordered model. The RG equations satisfied by the 
model, written in terms of the relevant thermal and magnetic fields t and h, 
with corresponding RG eigenvalues i/t and y/j, and the marginal dilution field 
ip, are given by ^ = y^r, ^ = ynh, ^ = q- gc,where b is the length 
rescaling factor and I = In 6. When q > Qc, the dilution field ip is relevant 
(and the phase transition is of first order), while in the regime q < Qc, ip 
is irrelevant and the system exhibits a second-order phase transition. The 
case q = qc is marginal. This picture is qualitatively correct, and in fact the 
critical value of the number of states which discriminates between the two 
regimes is qc = 4. In the q direction, qc — 4: appears as the end of a line of 
fixed points where logarithmic corrections are expected. At qc = 4, the RG 
equations were extended by Cardy, Nauenberg and Scalapino (CNS) [18, 19] 
and then by Salas and Sokal (SS) [20]. As a result of the coupling between 
the dilution field ip, and r and h, they were led to non-linear equations, 

dr 

= iVr + Vri^i^V, (16) 

= {Vh + yhi,'4')h, (17) 

dip 



d\nb 



gW- (18) 



The function g{ip) may be Taylor expanded, g{ip) — y^2ip (1 + ^^"0 + ...). 
Accounting for marginality of the dilution field, there is no linear term at 
q = 4. Comparing to the available results (for example the expression of the 
latent heat for <? > gc by Baxter [21] or the den Nijs and Pearson's conjectures 
for the RG eigenvalues for q < qc [10, 11]), the parameters were found to take 
the values = 3/(47r), y/i^ = l/(167r), ?/^2 = I/tt and y^a = -l/(27r2), 
while the relevant scaling dimensions are yr — v^^ = 3/2 and yu — 15/8. 

The fixed point is at r = = 0. Starting from initial conditions r, h, 
the relevant fields grow exponentially with I up to some r = 0(1), h = 0(1) 
outside the critical region. Notice also that the marginal field '0 remains of 
order of its initial value, V' ^ O{ipo). In zero magnetic field, imder a change of 
length scale, the singular part of the free energy density transforms according 
to 

/(V'o,r) = 6-^7(^,1). (19) 
Solving Eqs. (|16m8p leads to 

I ^ _1 Inr + ^ In (^G{^o. 0)) , (20) 

(for brevity we will denote v — 1/yr = |, /i = y^'y^ = \)- Note that G(V'o, "0) 

would take the value 1 in Ref. [19] and the value '^'^'^'^'^'f in Ref. [20]. We 

can thus deduce the following behaviour for the free energy density in zero 
magnetic field in terms of the thermal and dilution fields, 
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/(r>o)=r-ff^^^^^^^)"'/(l». (21) 

A similar expression would be obtained if the magnetic field h were also 
included. The other thermodynamic properties follow from derivatives with 
respect to the scaling fields. The quantity between parentheses is the only 
one where the log terms are hidden in the 4-state Potts model, and thus we 
may infer that not only the leading log terms, but all the log terms hidden in 
the dependence on the marginal dilution field disappear in the conveniently 
defined effective ratiofl Now we proceed by iterations of Eq. (|20p. and even- 
tually we get for the full correction to scaling variable the heavy expression 



V'o^, , ,^ / , , n 31n(-ln|r|)\ / 31n(-ln|r|) 

-G(V'o,V') = const X (-In |r|) 1 + ^ i / i ^ i 



V' .i_^^V 4 -ln|r| J\ 4 - In | 
CNS ^ ■ — ^ 



-1 



2 - In 



m 



SS 



4 (-ln|T|) J \ -hi|T| V-ln|r|2 



F(-ln|r|) 

where CNS and SS refer to the results previously obtained in the litera- 
ture [18-20] and F(— In |r|) is the only factor where non universality enters 
through the dilution field ^o- This allows to write down the behaviour of the 
magnetization for example 



M(r) = B_|r|i/i2(-ln|r|)-i/s 



^^31n(-ln|r| 



4 - In ItI 



4 -In r / V 4 - In r ' ^ ' ' 



-1/8 

(23) 



4 Numerical techniques 

In the Monte Carlo simulations we use the Wolff algorithm [22] for studying 
square lattices of linear size L (between L = 20 and L — 200) with peri- 
odic boundary conditions. Starting from an ordered state, we let the system 
equilibrate in 10^ steps measured by the number of flipped Wolff clusters. 
The averages are computed over 10^ — 10^ steps. The random numbers are 
produced by an exclusive-XOR combination of two shift-register generators 
with the taps (9689,471) and (4423,1393), which are known [23] to be safe 
for the Wolff algorithm. 

The order parameter of a microstate M(t) is evaluated during the simula- 
tions as 
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qN,JX-l , , 

. (24) 

g- 1 

where is the number of sites i with Sj = m at the time t of the simulation 
and m S [0, 1, (g— 1)] is the spin value of the majority of the sites. N = l? 
is the total number of spins. The thermal average is denoted M = (M). Thus, 
the longitudinal susceptibility in the low-temperature phase is measured by 
the fluctuation of the majority of the spins 

XL=/3((iV^)-(iV„)2) (25) 

and the transverse susceptibility is defined in the low-temperature phase as 
the fluctuations of the minority of the spins 

while in the high-temperature phase x+ is given by the fluctuations in all g 
states, 

X+ = ^E((^')- W)' (27) 

where A^^ is the number of sites with the spin in the state ^^. The internal 
energy density of a microstate is calculated as 

its ensemble average denoted as i? = (E) and the specific heat per spin 
measures the energy fluctuations, 

C = -/32||=/32((e2)-(E)2). (29) 

Our MC study of the critical amplitudes is supplemented by an anal- 
ysis of the high-temperature (HT) and low-temperature (LT) expansions 
for g = 4 recently calculated through remarkably high orders by Enting, 
Guttmann and coworkers [16,24]. In terms of these scries, wc can compute 
the effective critical amplitudes for the susceptibilities and the magnetiza- 
tion and extrapolate them by the current resummation techniques, namely 
simple Padc approximants (PA) and differential approximants (DA) properly 
biased with the exactly known critical temperatures and critical exponents. 
The LT expansions, expressed in terms of the variable z = exp(— /3), extend 
through z^'^ for the longitudinal susceptibility and through z^^ in the case of 
the transverse susceptibility. The magnetization and energy expansions ex- 
tend through z^^. The HT expansion is computed in terms of the variable 
V = {1 — z)/{l + [q— l)z). The susceptibility expansion has been computed 
up to v'^^ and the energy expansion up to v'^^. 
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5 Analysis of the magnetization behaviour 

For the sake of simplification of the notations, we group aU the terms con- 
taining logs in Eq. into a single function H{—\n\T\) = E{—ln\T\) x 
F(— In |r|) where 

Ei-ln\T\) 



3i„(-i„H)y/ +1 1 



4 -Inlrl / V 4-ln|r 



(30) 



The function E contains all leading logarithms with universal coefficients, it 
is known exactly while the function F needs to be fitted. We thus obtain 
a closed expression for the dominant logarithmic corrections which is more 
suitable than previously proposed forms to describe an observable (Obs.) in 
the temperature range accessible in a numerical study: 

Obs.(T) ~ Ampl. X \t\' X iJ«(- In |r|) x (1 + Corr. terms), (31) 
Corr. terms = a|Tp/^ + 6± kl + . . . , (32) 

where • and ft are exponents which depend on the observable considered, and 
take the values 1/12 and —1/8 respectively in the case of the magnetization. 
Here we stress that the inclusion of a correction in Irp/"^ seems to be necessary 
according to previous work of Joyce on the Baxter- Wu model [25, 26] (of 
4-state Potts model universality class), where the magnetization is shown 
to obey an expression of the form M(r) = i3_|r|^/^^(l -I- const x \t\'^^^ + 
const x 

Irl'i/^). The exponent 2/3 comes out from the conformal scaling 
dimensions of Dotsenko and Fateev [9] , and its presence is needed in order to 
account for the numerical results (see also Ref. [27]). Caselle et al. [17] also 
considered |r|^/'^ term to fit the magnetization. Here we also allow inclusion 
of a linear correction in 6±|r| to account for possible non linearities of the 
relevant scaling fields [1]. The next term in |t|''/^ will be forgotten. 

In Fig. [3] we plot effective magnetization amplitudes B^fflr) vs |t|^/'^. 
From the available data for the magnetization, we define the following quan- 
tities, 

BcNs{r) = A/(T)|T|-Vi2(_in|r|)V8, (33) 
Bssir) = A/(r)|r|- 1, |,|)i/8 (^1 _ ^^^^^^^) " , (34) 
Be{t) = Af (r)|r|-i/i2£;i/8(_ jj, |^|)^ (35) 
which are expected to behave according to the corrections to scaling 

B_{l + a\T\'''^ + b±\T\ + ...) (36) 
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Fig. 3. Effective amplitudes (from MC data) as deduced from different assumptions 
for the logarithmic corrections. Symbols correspond to MC data, dashed and dotted- 
dashed lines are fits as explained in the text. 



The numerical results are shown in Fig. [3l From bottom to top the various 
symbols indicate the effective amplitudes with "no- log" at all, then with 
the CNS and the SS corrections and finally the effective amplitude where 
the known universal logarithmic terms have been included. The dashed lines 
correspond to a rough determination of the corrections to scaling including 
only the terms in alrp/"^ in the limit |t| 0, and the dot-dashed lines 
include also the terms in 5|t|. From this plot, we deduce that none of the 
three effective amplitudes in Eqs. \33f^35^) can be correctly fitted by Eq. ^30fl . 
since the coefficients of the correction terms (e.g. the coefficient a which is 
estimated directly by the slope at small |t| values) strongly depend on the 
range of fit. This undesirable dependence of the coefficients on the width of 
the temperature window is shown in the first six lines of table [TJ 

In order to improve the quality of the fits, one has to take into ac- 
count the correction function i^(— In |r|) and to extract an effective function 
Ff,ff{— In |t|) which mimics the real one in the convenient temperature range. 
This is done by fitting B^f / (t) to a more complicated expression, 

which means that we include the corrections to scaling and the non universal 
function function F{— In |r|) taking the approximate expression 



/ , , ,^ C-i C2ln(-ln|T|) 



(38) 
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Table 1. Fits of the effective amplitude of the magnetization. 



Beffir) 


|rj^'' ■'-window 




a 


b 


BcNs{r) 


[0,0.15] 
[0, 0.45] 


1.07 
1.05 


-0.98 
-0.66 


0.94 
0.29 


Bss{r) 


[0,0.15] 
[0, 0.45] 


1.11 
1.10 


-0.77 
-0.45 


0.56 
-0.06 


BE{r) 


[0,0.15] 
[0, 0.45] 


1.14 
1.13 


-0.47 
-0.25 


0.16 
-0.27 


Bef{t) 


[0,0.15] 
[0, 0.45] 


1.16 
1.16 


-0.20 
-0.18 


0.02 
-0.02 



While a and h are coefficients of corrections to scaling due to irrelevant opera- 
tors, Ci and C2 are effective coefficients of logarithmic terms which, in a given 
temperature range, mimic a slowly convergent series of logarithmic terms de- 
pending on a non universal dilution field. Therefore, we expect that different 
fits made in different temperature windows will produce different values of 
Ci and C2 while a and b (and of course also the magnetization amplitude 
should be relatively less influenced by the window range. The choice of 
values for Ci and C2 is thus partially arbitrary and the values quoted should 
be specified together with the temperature window where they are appropri- 
ate. In the following, we obtain Ci ~ —0.76 and C2 — —0.52 in the window 
|-r|2/3 g j-Q^ 0.35], which yields an amplitude i?_ ~ 1.157. The resulting a and 
b coefficients now appear very stable. This is checked in Fig. |4] where the 
quantity 

BEF{r)^M{T)\T\-^'^^[E{^\n\T\)F{-\n\T\)f'^ (39) 

is reported toghether with the previous curves and fitted as indicated in 
table [T] As an independent test, we add the SE data which are superimposed 
to the MC data at small values of jrj only for this latter assumption of 
effective amplitude. 

Eventually, our approach confirms expression Eq. ((23|) for the magnetiza- 
tion, with the function F(— InjTj) given in Eq. ([55)1 and the parameters Ci 
and C2 given above for the appropriate temperature window. Nevertheless, 
we have to stress that the different effective amplitudes should all reach the 
same amplitude B- in the limit jrj 0, since Bcns{t) is in fact an approx- 
imation of Bss{t), which is an approximation of Be{t), which eventually 
approximates Bef{t). An attempt of illustration of this behaviour is shown 
in Fig. [5] where dotted lines (which are only guides for the eyes) all converge 
towards the unique value i?_ ~ 1.157. Here we stress that we have deleted 
the points of SE data which are too close to the critical point, since the series 
are no longer reliable because the current extrapolation procedures are in 
principle unable to approximate the complicated structure of the singularity 
involving log corrections. 
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Fig. 4. EfTective amplitudes (Irom MC and SE data) as deduced Irom different as- 
sumptions for the logarithmic corrections. The upper curves correspond to Bef{t) 
and the thick solid lines to SE data. 
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Fig. 5. Zoom of the effective amplitudes in the vicinity of the critical point. MC 
data (symbols) and SE data (thick solid lines). 
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6 Universal combinations for the 4-state Potts model 
and conclusions 

The other quantities can be analyzed along the same lines. The important 
point is that now the fmiction F^f / (r) is fixed in the corresponding tempera- 
ture scale and thus the remaining freedom for the other physical quantities in 
Eq. ((3T|l is only through the leading amplitude and the coefficients of \t\'^^^ 
and |r|-terms in the corrections plus possibly the background terms. It is still 
a complicated task to perform this analysis, but the current results for the 
universal combinations mentioned in the introduction appear in the following 
table. 

Table 2. Rough estimate of the universal combinations of the critical amplitudes 
in the 4-state Potts model. 



A+/A. 


r+/rL 




Rc 


source 


1.° 


4.013 


0.129 


0.0204 


[12,14] 




3.14(70) 




0.021(5) 


[17] 




3.5(4) 


0.11(4) 




[16] 


1.00(1) 


6.7(4) 


0.161(3) 


0.0307(2) 


here 



" exact result from duality 

Our work is "one more" contribution to the study of this problem and 
brings some answers, but also raises new questions. Indeed, our results dis- 
agree with previous estimates, but we cannot claim for sure that our estimates 
are more reliable than those of other authors. What is extremely clear is that 
the groups who studied numerically universal combinations of amplitudes in 
the 4— state Potts model all noticed the extreme difficulty to take into account 
properly the logarithmic terms. We believe that our protocol is self-consistent 
in the sense that our criterion is to obtain a relative stability of the correc- 
tion to scaling coefficients. The results that we report here, although a rough 
estimate which calls for deeper analysis, reach a reasonable confidence level. 
If this is indeed the case, one should identify the reason of the discrepancy 
from the theoretical predictions of Cardy and Delfino. In the conclusion, and 
in a footnote of one of their papers, Delfino et al [14] (p. 533) explain that 
their results are sensitive to the relative normalization of the order and dis- 
order operator form factors which could be the origin of some troubles for 
the ratios Fj^/Fl and Rc- This possible explanation seems nevertheless to 
be ruled out (as mentioned by Enting and Guttmann already) by the very 
good agreement between the theoretical predictions and all numerical studies 
(both MC and SE) in the case of the 3-state Potts model. Eventually let us 
mention that the two-kink approximation used by Cardy and Delfino is exact 
for q = 2 (Ising model) and quite good for q = 3, but probably questionable 
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close to the marginal case g ^ 4. As a conclusion, we are afraid that this 
work opens more questions than it bring answers. 
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